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A. Notzhon 


Oe AY} mesas: soinesd co by Accumutato: (low! end Y resizes 


gh). 


< (A~Z.F.=$X} means: Accumulator jas to ccotain, cr contains the 
cuntents of location $X. IE $X equals zero, the zero t'ag iZ.F.) is set or 
Tuc3i be set, otcerwise che zero tag is or muse be clear. 
L. Remarks 
1. “or some coutines presented below, the entry anc/or exit values o’ 
"he accumulator, the X register and the Y register cre given. u nc 
entry it specified, no entry is necessary. L no exit or no: ai! reg sess 
=f an exit are specified, the registers ncs specified may have uay.<- 
eilte v4 vels7s aiter the exece.ioa of tae Applesoft rourne. 
& Fc. each reuaine, the memo.y locations that may be modifies by the 
feror free) execution, of the rovt nea given. 
“` Waraing 
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When work’ nyg with m.l programs Mat are called from BASIC, one 
may wish to use zero page locaaons fa stor temporary resolts. 
However, a nuisbec ol zero pase locations gie mt2zelized to cettam 
values @ the roli/ wann stan of Ap plesott and ch n gi: ihe vans 
Tents af chee ba anens may lead 10 unexpected resultè byrcher 
More there aw d number ot locations in which ariei: sTures 
idormation during the eaceation ca the BAS prograta, suck asthe 
cunnt line number or the pointer co che hne from whick data is 
doing Mad Clovberog one of these locations usuaily has ¿e ettect 
hat the progran will ciash sooner or later 


In ame toe svol problems wher. working with feru nape 
oddiesses itas therefore recommended (se consult the zero paye uage 
Ma» ain tbe Applesoit manuail tise (See also bie taemoury atias con: 
tucca by Viol W F. Lusobert, published in the August 1979 issue 


of MICRO} 
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At that time, however, i discovered 
that the Appie can also be programmed 
rather easily im machine language with 
the heip of the mint-assembler Since | 
was interested to know what speed gain 
could be obtained, 1 translated the per- 
mutation program in machine code To 
uty surprise, the programs tan about 675 
times faster (approximately 13,500 per- 
mutations per secondi! than the BASIC 
permutation program. Of course, | 
immediately got my pentonnno pro- 
gram and translated this m machine 
code too. The 2,339 salutions now 
cane outan less than 3 bouts, which 
acar also a considerable gainan speed 
as compared to the BASIC program. 


When to Use Machine Language 
Veograms or Subroutines 


Some programs, like those men- 
ioned above, can vasily be translated 
~om 3ASIC to machine code. 
However, fcr che majority of the pro- 
e ms that I intend to write, this 1s not 
the o.:2, since in these programs 
floating point variables rather than 
"one byte” variables have to be used. 
For some floating point arithmetic such 
as sddition and multiplication, it is 
probably possible to write the routines 
yoursel®, but for functions such as the 
sine and the loganthm this would 
mean alot of work. Furthermore, being 
busy with trying to reinvent the 
wheel''is sot a very stimulating idea. 


However, there is a fairly easy way 
out of this problem. All the routines 
needed for floating point arithmetic, 
have to be somewhere in the Appiesolt 
ACM, so ail one has to do is list 
oo,lesoft and try to understand how it 
verns. Afar locating the entnes of the 

Le VP pout reuanes. these routines 
on oan be called ce the machine 
vanguare md} program. Although the 
de ssocess con be writter down ina 
IEW NES, - wet Seven Werks of 
hard work vers b knew enouyh ot 
Applescts to write, as an exe cine, R 
nama mulnpheanon subroutine m 
wl whieh can be called trom BASIC 
ty means of the & sembol This 
program runs about 8 times taster than 
a BASIC matrix multiplication 
subroutine and further has the advan- 
tage that the names of the matrices can 
be passed in an easy way. On the other 
hand, a disadvantage is that the m.i. 
program uses more memory space than 
the BASIC program. At the end ef this 
article, the matrix multiplication 
program will be more extensively 
discusse?. 


a) 
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Tabie 1: Applesolt Routines (continusd) 


Of course, this warning does not apply to {most ot the! zero page 
Jocanions that may be modified by the routines described be.ow. For 
instance, if one uses neither the power fuarction nor SQR nor 
triyonoractical tuncions, at will be sate to ure lucations $SA-$H8E, 
since these locations are used by none of the other functions 
{routines} listed in this table. 


if neither strings nor high-resolution graphics nor ON ERR 
statements are used, one can (probably) salely store temporary 
results in the following zero page locations. 


$6 $9, $17-$1F, $58-$5D, $71-$72, S$CE-$D5, $D7, $D9-SEF, 
$P4 SFF 


I. Description and entries of the routines 


A. Charget-Charcheck 


l. 
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Purpose 


The memory locations $B8 ana <B9 coatain—during the execution 
of « BASIC program—a text pointer waich points to the lest re- 
trieved character of the BASIC program. The Charget routings can be 
used to load the next character or the current character again’ in the 
Accumulator. To determine whether the character equals a predeter- 
mined syinbol one of the Charcheck routines may b> used. 


Charget routines 


$Bl: Advance text pointer and load next character in the 
Accumulator {spaces are ignored]. 


Exit(A = next character,X = entry, Y = entry}. 


Exit Status: Carry is clear if character is a digit (hex 
value:30-39), otherwise carry is set. Zero tlag is set if character 
equals 0 { = end of line sign} or 3A | = end of statement sign, 1.¢. 
1) otherwise zero flag is cleared. 


Modifies $B8,$B9. 

¿B?: Load current character another time in the Accumulator. 
Exit(A = current character,X = encry, ¥ = entry). 
For status see subroutine $B1. 


Charcheck routines 
$£97D: Check whether character in Accumulator is a letter. 
Entry{A],Exit(A «ertry,X =entry,Y = entry}. 


Brit Status: Carry is set if character 1s a letter, otherwise 
carry is cleared. 


1 e following 4 routines can be used to check whethe: the text 
pornter pnts to a speciiic symbol it the resule ot the cnech ts 
powtive, the next character as loaded an the Accumulator by means 
ol the execution of subroutine 381 In the other case, the messaec 
MOVNTTAX .RROR” as displayed and Applesolt returns to BASIC 
command level. The exits of the 4 routines are: 


Exits(A = next character,X = entry, Y =O), modify $B8,$B9. 


$DECO: Check whether the character that is pointed to by the text 
pointer equals the character in the Accumulator. € e 


Entry jA). 


$DEB8: Check whether the text pointer points to a ngnt 
parenthesis. 


SDEBB: Check whether phe text pointer points to a left parenthesis. 
$DEBE: Check whether the text pointer points to a comma. 
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B. Compare 


1. Purpose 


The compare routines can be used for comparing a real variable in 


the MIY 


with areal variable in the SEP or a real variate in memory. 


2. Compare routines 


SLFG6A: 


Type of Comparison 


ev 


AAN 
v^ 


$EBB2: 


C. Conversion 


1. Purpose 


Compare MIP with SEP according to the status of the com- 
parison in location $16. The result of the comparison (1 if 
truc, O if false} is converted to a real variable in the MEP, 
The various types of comparisons are listed below. 


$16 has ta be put 
equal to: 


Result comparison 


if SEP > MEP, else 0 
if SFP = MFP, else 0 
if 


if SFP = MFP, else 0 
1 if SFP # MFP, else 0 
if SFP s MFP, else O 


AWG & Vv = 
Ds jee pun hd t pa 


Modifies $60,$61, MFP, SFP. 
Compare MFP with memory [> A, Y}. 


Entry(A,Y}, Exit(A=FF if MFP < memory, A«0 if MFP = 
memory, A=1 if MFP > memory], modifies $60,$61. 


The Conversion routines can be used to convert: 


a} areal 


in the MFP to an integer 


b) a one or two byte integer to a real in the MFP 


Unless specified otherwise, all integers are assumed to be two's 
complement integers. 


2. Real to integer conversion routines 


S$EBF2: 


$752: 


SE10C: 


$E108: 


$DA65: 
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Convert MFP to integer. T..c number in the MFP must be 
between -23' and 20 /notation: - 29! < MFP < 233). 
Result is stored in mantissa ot MFP (locations $9E-$A1). 


Exit(Y = 3}, modifies MFP. 


Convert MFP, where —2'* < MFP < 2'4, to two byte 
integer. Store result in $50 {low} and $51 (high]. 


wee fA 200, Y a yG, Modes $90,$51 MEP 


Remark. "Wrap around” occurs if the absolute value of the 
number inthe wth? ts larger chen 2'* ~ L. 


Convert MFP, where ~2'§ < MFP < 2, to wwo byte 
integer. Store result in $A (high) and $A1 {low}. 


Exit (Y= 0), modifies $60,$61,MFP. 


Same as $E10C, except that entry-value of MFP must be: 
0 =< MFP < 2'5, 


Pack extension byte in MFP and convert MFP, where -2'3 
< MFP < 2'5, to two byte integer. Store integer (high 
byte first} in (> $85,586). 


Exit|Y = 1), modifies $60,$61,MEP. 


(continued) 


MICRO - The 6502 Journal 


=e ey 1 maie- pee nan eaten) 





ere SSS rs psp ich SSP ae A SS 


aaa aT 


An important point to note is that, 
as a conscquence of using floating print 
arithmetic, there is a signiticant drop of 
the speed gain, namely trom a factor 
675 sbtained with the permutation pro- 
Kram to a factor & obtained with the 
Matrix miuluply atan program The 
treason is that- when mulcplving 
Matices -a relatively large pornon of 
the CPU time ais used for the 
multiplication and addition of floating 
point numbers. Whether this is done 
under control of a BASIC program, ot 
by calling the appropriate roudnes in 
Appicsofe from a m i. program, makes 
no difference, since in both cases the 
same multiphcavon and addiuon 
routines are used, The gain of speed 
that occurs in the m.l. matrix 
multipiication program is obtained by 
spo t-cutting the time-consuming 
determinatica of the pointers to array 
elements in BASIC. t. 

It will now also be clear that it does 
not make any sense to caiculate for ia- 
stance, 1000 logarithms by means of a 
m.l. program. When written in BASIC, 
thus 


10 FOR! = 1 TO 1000: 
A = LOG (|): NEXT 


the program will run approximately 23 
seconds. About 90% of this time, the 
computer will be busy with the calcula- 
tion of the logarithms, and about 10% 
of the time with the parsing of the 
statements and the evaluation of the 
FOR...NEXT loop. When writing a m.i. 
program to calculate the logarithms, 
one may expect it to run no more than 
10% faster than the BASIC program, 
since as to the calculation of the 
logarithms, no time can be saved. 


Therefore, with respect to gaining 
speed, it is only profitable to write a 
m.l. program or subroutine if, in this 
way, time-consuming access to array 
elements can be shortcucted or 
iHerdave parts of the program: can be 
made more efficient: Some examples 
where m l. routines will be useful are: 
finding the largest element of an array, 
calculating the inverse of a matrix, 
sorting the elements of a vector, or 
calculating probabilities under a 
bivariate (log) norma! distribution, 


A parc from gaining speed, there may 
how ever be other arguments fer writing 
ro.l. routines. For instance, one may 
wish to extend tzpe or disk versions of 
Appiesoft with some self-written 
BASIC commands or functions. Also, it 
can be attractive to make frequently 
used subroutines more independent of 
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Table 1: Appiesoft Routines (continued) the main program, so that parameters 
can be passed by value rather than by 
3. Integer to real conversion routines name, which in BASIC is only possidle 
wal ; ; , ; by means of a lot of PEEKs and POKEs. 
$E2F2: Convert two byte integer in A (high) and Y (low} xo real in Last but not Icast, one may like the 
MFP. challenge involved in wnting m l. 
Entry{A,Y}, Exit(Y <0], modifies MFP, puts $11 equal to programs. 

zero. . 

The Main and Secondary Floating 


eae i ae i | 
$E301: Convert one byte integer in Y to positive real in MFP. {The Point Accumulator | 
t 
I 


integer in Y is thus not interpreted as a two's complement 


integer. | Before presenting the Applesoft 


Entry{Y}, Exit(Y =0), modifies MFP, puts $1) equal to zero. routines that can be of help when 
writing ml peogeams, the main and 


$EB93: Convert one byte sateger in Accumulator to real in MEP. 
secondary floating port accumulator, 


Entry(Aj, Foat(Y 0), modifies MEP. {heneetorth to be abbreviated as MEP } 
$DEES. Pull anteyer (%) variable from memory (®$A0,$SAL1) into A and SEP respectvely], will shortly be | 
jhigh} and Y (low). Next, convert integer to real in MEP. H OS a the arithmenĘsal 
) Sa and mathematical routines use the { 
Exit(Y = 0), modifics MFP, puts $12 equal to zero. MEP and/or the SEP. The MEP oc- 
ae cupies she memory locations $9D-$A2 } 
-p7 and $AC. The exponent of the floating 
point number is in SYD jin excess 80 
-- Purpose code], the mantissa is in $9E-$Al, and i 
= Tae Copy routines can be used to its sign is in SA2. Location SAC is used : 
a) pull data (from memory) into the MFP or the SFP in most floating point routines as an 
b) pack the MFP and store the MFP in memory extra mantissa byte, to increase the 
c] copy the MFP into the SFP and vice versa presen of the saleulavions, This loca- 
d) push the MFP on stack or pull the SFP from stack ae a eA eng the aaa 
l : x (e) ow one ca 
The Copy routines are for real variables only. For routines that ailal d the quate of the MEP to a 
handle integer (%} variables see Conversion. decimal number is given in example 1 
; (page 31). The sign of the number is 
2. MFP routines positive, since the first bit of $A2 is 
SEAF9: pull memory (- A,Y] into the MFP and put the extension zero. In case this bit equals one, t: 
byte equal to zero. sign of neal in the Sheik e 
; ` Z vgs negative. The exponent is calculated by 
Pen eae aca ~$9D,X «entry, ¥=0}, modifies converting the hex number öl in $9D 
rare Á to decimal, which gives 129, and by 
$EAFD: Pull memory (-» $5E,$5F] into the MFP and put the exten- subtracting the excess |=80 {bex} or 
sion byte equal to zere 128 {decimal} } from it. The method 
Exit[A = Z.F. = $9D,X =entry,Y =O], modifies MFP. that is used to convert the mantissa to 
ae an i decimal is essentially the same as the 
$DE10: Pack extension byte in MFP and push MFP on stack (6 method used to convert a normal hex 
bytes}. number to decimal, except that instead 
ExitjA=Z.F.=$9D}, modifies $5E,$5F,MFP. of E onsen S m 256, 4,096, 
Ld LA ‘ 
= The following fcur routines pack the sign and the extension byte Ste. the reciprocals: of these numbers ' 


in the MFP, store the MFP in the locations indicated and put the 


have to be used. | 
extension byte equal to zero. 


The number zero forms an excep- 
tion to the miles mentioned above. 


For. ‘cur routines the exits are: Applesofe considers a number to be cero 


ExitysAeZ E $9, Y= 0), modify $5E,$5F, MFP. if the exponent ($9) equais zero, an- 
$EB1E: store MEP in $58 $9C: dependent of the value of the mantissa. 
$EB21: store MFP in $93-$97 The cesults trom arithmetical 
$EB27: store MFP in (> $85,86) operations and matheraatical functions 


in Applesoft are, in general, placed in 
the MFP. Next, the MFP is usually nor- 
; malized and pushed on the stack or 
3. SFP routines stored in memory. The normalizing of 
$E9E3: Pull memory {> A,Y) in the SFP and determine $AB { ~ the the MFP means that the bytes or the 
exclusive OR of the signs of the numbers in the MFP and mantissa are rotated to the left (zeros 


$EB2B: store MFP in (> X,Y} 


the SFP). enter at the right) until the left-most 
Entry[A, Y}, Exit (A=Z.F. =$9D,X =entry,Y =0], modifies | bit of $9E equals one. At every rotation 
$5E,5F,SFP,SAB. the ie is decreased W7 one, since 

a rotating the manti itt 
$EIE?: Pull memory (-> $5E,$5F} in the SFP and determine SAB. li Ea E SU O ORe 


means multiplying the number in the 
E it[A = Z.Fr= $9D,X »entry,Y = 0), modifies SFP,$AB. MFP by two, and this number must, of 
' (continiied) course, remain the same. f 
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Uf, after the normalizing process, the 
MFP has to be stored in memory, it 
must be packed because the MFP cc- 
cupies 7 bytes of memory, whereas 
Applesoft reserves only 5 bytes for the 
storage of real variables. In the packing 
routine, first the mantissa 1s rounded 
off by considering tho left-most bit of 
the extension byte. If this bit equales 
one, the mantissa is incrcased by one, 
otherwise the mantissa remains the 
same. Then the sign is packed into the 
floating point number. If the sign is 
positive, the left-most bit of $9E is put 
equal to zero, otherwise it remains 
equal to one. Note that the sign can be 
packrd in this way because the first bit 
of $9E contains no information since it 
always equa's one alter normalizing. 


The SFP occupies the memory loca- 
tions $A5-SAA. Tne exponent is in 


$A5, the mantissa in $A6-SA9, and its ` 


sign in SAA. The SFP has no extension 
byte. For the arithmetical and 
mathematical operations requiring two 
operands, the first operand has to be 
put in the MFP and the second operand 
in the SFP. Thus, loading the SFP and 
the MFP with two numbers and doing a 
JSR to, for instance, the multiplication 
routine, leaves the product of the 
numbers in the MFP. For some 
arithmetical routines it is necessary to 
determine—before the routine is 
executed—the exclusive OR of the 
signs of the numbers in the MFP and 
the SFP. The result must be stored in 
location $AB. This implies that the 
` first bit of SAB must be one if the signs 
differ, otherwise the first bit has to 
equal zero. However, in most cases the 
user does not have to bother about 
determining the value of $AB, since it 
usually is not necessary to load the 
MFP and/or the SFP "by band.” 
«pplesoft provides us with a lot of 
routines that can be used to get floating 
point numbers from memory, unpack 
“Bort, 2. i alace thers ic ths MFP or che 
SFP. All the rouanes that puli memory 
in the SFP ilsu set $AB to the right 
value, 


The Use of Applesoft Routines 


The Applesoft subroutines that are, 
in my opinion, the most useful for m.l. 
programmers are listed in table 1. A 
distinction has been made between 
various types of subroutines, such as 
Copy, Errors, Conversion and 
Mathematical routines, etc. Rather 
than discussing each of the routines 
separately, a (very) simple examoie 
will be given to illustrate how to work 
witb them. For a good understanding of 
this example, it ir advisable to read the 
„eneral remarks ic table 1 first. 
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Table 1: Applesoft Routines (continued) 


$DE47: Pull stack in the SFP and determine $AB. This routine will 
usually be used in combination with subroutine $DE10. Ia 
that case it is for a successful execution of routine $DB47 
necessary to push the retum address of SDE47 on stack 
{high order byte first) before executing $DELO. Contrary to 
most other routines described here, §$O847 mus: be ex: 
ecuced by means of e JMP instruction. 


Exit(A=Z.P.=$9D,Xeentry,Yeentry), modifies 
SPP,SAB. 

SFP/MFP routines 

$EB53: Copy SFP into MFP, put extension byte equal to zero. 
Exit(A=$9D,X =0,¥ sentry}, modifies MEP. 


SEB63: Pack extension byte in MFP and copy MFP into SFP, put ex- 
tension byte equal to zero. 


Exit(A = $9D,X =0,Y sentry}, modifies MEP, SEP. 


SEB66: Copy MFP (withou: sitensica byte) into SFP, put exten- 
sion byte equal to zero. 


Exit(A ~$9D,X =0,Y = entry), modifies MFP,SFP. 


E. Errors 


1. 


Purpose 


If an error is detected in a m.l. program, one of the error zoutines 
may be used to print an error message. 


. Error messages 


To print an error message, load the X register with the code of the 
message and execute a JMP to $D412 or execute a JMP to one of the 
locations listed behind the error messages. After printing tbe error 
message, Applesoft returns to BASIC command level {unless an ON 
ERR statement has been executed). i 


Error message JMP location 
NEXT WITHOUT FOR $DDOB 
SYNTAX ERROR $DEC9 
RETURN WITHOUT GOSUB $D979 
OUT OF DATA ~ 
ILLEGAL QUANTITY $E199 
OVERFLOW SE8DS 
OUT OF MEMORY $0410 
UNDEF'D STATEMENT $D97C 
BAD SUBSCRIPT $E196 
REDIM’D ARRAY — 
DIVISION BY ZERO SEAE? 
ILLEGAL DIRECT $E30B 
TYPE MISMATCH $DD76 
STRING TOO LONG _ 
FORMULA TOO COMPLEX $£430 
CAN'T CONTINUE _ 
UNDEF’. FUNCTION $E30E 


F. Expressions 


fe 


ANCRO ~The 6502 Journal 


"ma eee 


Purpose 


The Expressions routines can be used to evaluate expressions ‘in an 
& statement. When calling an expression evaluation routine, the 
text pointer in $B8 and $B9 must point to the first character of the 
expression. After control is returned from the evaluation routine, the 
text pointer points to the first character behind the expression. In the 
evaluation routines below, this character is called the terminal sign. 
The terminal sign might, for instance, be a comma, but also a special 
character such as a "4#". The locations that are modified by the 
routines are not specified here, since these depend on the type of the 
expression. ` 


A 
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2. Expression evaluation routines 
$DD67: Evaluate expression to next terminal Sign, store result in 


$E105: 


SEGFS: 


G. Init 


MFP. 


Evaluate expression to next terminal sign and convert 
result, which must be non-negative, to two dyte integer in 
$A0 {high} and $A1 (low). 

Exit(Y = 0}. 

Evaluate expressiim to next terminal sigu and convert 


result, which must be non-negative, to a one byte integer in 
SAI. 


ExitiA = terminal sign,X = $A1,Y =}. 


1. Purpose 
initialize mantissa of the MEP or the SEP. 


2. Initialization routines 


$EC40: Init mantissa MFP {except extension byte] and Y to value in 


SE84E: 


Accumulator 

Entry(A), Exit(A =entry,X = entry, Y = A}, modifies MFP. 
Put MFP {$A2 and $9D} equal to zero. 

Exit(A = Z.F. =0,X sentry, Y= entry), modifies MFP. 


H. Mathematical ! (rputines with one operand) 


SEBAF: 


SFOSE: 


SEEDO: 


SEFEA: 


SEFO9: 


$DE98: 


SEB9O: 


$EB82: 


SEFF): 
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MFP = ABS{MFP} 
Exit(A »entry,X = entry, Y = entry}, modifies MFP. 

MFP = ATN{MFP} i 

Poe SSE, $S5F,$62-$66,$92-$9C, MFP,$A3,SFP,SAB,SAD, 


MEP = -MFP . 
Exit(X =entry,Y = entry}, modifies MFP. 
MFP = COS(MFP} 


Modifies $D,$16,$5E,$5F,$62-$66, $92-$9C, MFP,$A3,SFP, 
$AB,SAD,SAE. 


MFP = EXP(MFP} 


Modifies SD, SSE, $5F,S62-$65,$92,$98-$9C, MEP, $A3,SFP, 
$AB,SAD,SAE. 


: MFP = INT(MEP} 


Modifies $D,MFP. 
MFP = LCG(MFP} 


Medifies SSE, S5F,562-$66,$92-$9C, MEP, $A3,SFP,$AB,SAD, 
SAE. 


MFP = NOT({MFP}. This routine returns MEP = l if MFP = 0, 
else routine returns MFE = 0. 


Modifies MFP, puts $11 equal to zero. 
MFP = SGNiMEFP} 

Exit{Y =0}, modifies MFP. 
Accumulator = SGN(MFP} 


-Exit(A=FF if MFP <0, A=0 if MEP = O and A=! if 


MFP > 0,X sentry, Y = entry). 
MFP = SIN(MFP} 


Modifies §D,$16,$5E,$5F, $62-$66, $92-$9C,MFP,$A3, SFP, 
$AB,$AD,$AE. 
(continued) 
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Suppose one wishes to translate a 
BASIC subroutine to m.l. In that case 
the m.l. routine can be called from 
BASIC by means of the & symbol, The 
& symbol causes an unconditional 
jump to location $3F5 where the user 
can insert a IMP instruction to the start 
of the m.l. program. í 

e 

After the execution of the & symbol, 
the text pointer af BASIC, which Is in 
the locations $B8 and $89, points to 
the next character of the line (spaces 
are ignored). Thus, if we have :he line 


10 & A1, BQ. C 


where Ai, BQ and C are reals, she text 
pointer points—after the execution of 
the & symbol—to the A. Suppose we 
wisk «o multiply Al and BQ and store 
the zesult in C. We then first have to 
determine the starting locatior of the 
storage area of the value of Al in 
memory. This can be done by making 
use of the subroutine $DFE3, listed 
under the heading Names in table 1. A 
JSR to SDFE3 in the m.. program ex- 
ecutes an Applesoft routine which puts 
the name of the variable jin tkis case 
Al) in $81 and $82; the status of the 
variable {in this case real) in $11 and 
$12; the pointer to the location of the 
variable in $9B and $9C; and, most im- 
portant, the pointer to the value of the 
variable in $83 and $84, as well as in 
the Accumulator and the Y register. 


Now that the starting location of 
the value of Al is known, the value of 
Al can be pulled into the MFP. For this 
purpose, the Copy routine SEAF9 can 
be used. Since the entry of this routine 
corresponds with the exit of SDFE3, the 
subroutine call to SEAF9 can be placed 
directly behind the subroutine call to 
$DFE3. 


Now that we have stored Al in the 
MEP, we can proceed to analyzing line 
10. After the execution of subroutine 
SDFE3, the text pointer points to the 
first character behind the name of the 
variable, which is—in our example—a 
comma. lf one plans to write a serious 
m.l. program, it might be useful to 
check whether there is indeed a comma 
behind the name. 


For checking purposes, various 
routines are listed under the heading 
Charget-Charcheck. For instance, to 
check whether a comma is present, a 
JSR to SDEBE can be executed. In case 
the character is not a comma, the 
“SYNTAX ERROR" message is 
displayed and Applesoft gives a warm 
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start on BASIC. If, on the other hand, a 
comma is present, the text pointer is 
advanced an! points now to the letter B. 


To obtain the starting location of 
the storage area of the value of BQ, 
again a JSR SDFE3 is executed. Since 
Al and BQ have to be multiplied, BQ 
must be stored in the SFP To ac- 
complish this, subroutine $FYE3 is 
used, which also can be placed directly 
behind the JSR $SDFE3 instruction, 
because the exit of SIOFE3 corresponds 
with the entry of SE9E3 Note that itas 
necessary to fill the MFP before the 
SEP, because SAH is set when HQ is 
pulled in the SFP. 


As cap Se seen, the entry of the 
toultiplication routine $E982 cor- 
responds wth the exit of $E9E3. Sc 
after the JSR to SEVE3, the multiplica- 
tion can be carried out by means of a 


_. ISR $E982. Note that the m.|. program 


~~ 


can be reduced by several bytes by 
using JSR $£97F instead of the last two 
mentioned subroutine calls. 


Finally, the result of the multiplica- 
tion, which is in the MFP, has to be 
stored in C. Before this is done, a JSR to 
$DEBE is executed to check whether 
the text pointer points to a comma. 
Next, the starting location of the 
storage arca of the value of C is deter- 
coined by means of a JSR SDFE3 in- 
structior.. To store the value of C in 
memory, the Copy routine $£B2B can 
be used. Since the entry of this routine 
is X,Y), whereas the exit of $DFE3 is 
(A, Y], the instruction TAX must be in- 
serted before the instruction {SR 
$EB2B. 


Àf er the last execution of SDFE3, 
cae text pointer points co the end of line 
10, soa RTS instruction in the m.l. pro- 
‘Tair turns contre! io the BASIC pro- 
gum which will restart execution at 
the Lane number following line 10 The 
complete m l. program is given in 
example 2 (page 40}. 


The routine $DFE3 can also be used 
to find the start of the storage area of 
integer (%} variables, elements of 
arrays, and arrays. If one wishes to use 
matrix expressions in the & statement, 
it is necessary to store the hex value 40 

in $14 because otherwise Applesote 
will interpret the matrix names in the 
& statement as names of simple 
variables. Be sure you don’t forget to 
put $14 back on zero before returning 
1o BASIC, because otherwise strange 
‘hings may happen. 


I. 





Table 1: Applesoft Routines (continued) 


$EE8D: MFP = SQR(MFP} 


$PO3A: 


$EFAE: 


Modifies $D,$5E,$5F,$62-366,$8A-$8E,$92-$9C,MFP,SA3, 
SFP,SAB,S$AD,SAE. 


MEP = TAN{MFP] 


Modifies $0,$16,$SE,$SF,$62-$66,$84-$8E,$92-$9C,MEP, 
$A3,SFP,$AB,$AD,$AE. 


MFP = RND(MFP). See Applesoft manual for argument RND 
function. 


Modifies $5E,$5F,$62-$65,$92,MFP,SFP,$C9-$CD. 


Mathematical It (routines with two operands} 


Add 

$E7C1: MFP = SFP + MFP, $AB must be determined before 
subroutine call. 
Entry[A = Z.F. = $9D), mcdiries $92,MFP,SFP. 

$E7BE: Pull memory j> A,Y} in SFP, determine SAB, add: MF? = SFP 
+ MFP. 
Entry(A,Y}], modifies $5E,$5F,$92,MFP,SFP,SAB. 

AND 

$DF55: MFP = SFP AND MFP. Routine returns MFP = 1 if MFP and 
SFP are both unequal to zero, else routine retums MFP w O. 
Modifies MFP, puts $11 equal to zero. 

Divide 

$EA69: MFP = SFP/MFP, $AB must be determined before subroutine 
call. 
Entry{A=Z.P. = $9D), modifies $62-$66,MFP,SFP. 

$EA66: Pull memory (-» A,Y} in SFP, determine $AB, divide: MFP = 
SFP/MFP.. 
Entry(A,Y), modifies $5E,$5F,$62-$66,MFP,SFP,SAB. 

Multiply 

$E982: MFP = SFP x MFP, $AB must be determined before subroutine 
call. 

_ Entry{A=Z.F.=$9D), modifies $62-$65,MFP. 

$E97F: Pull memory (> A,Y} in SFP, determine $AB, multiply: MFP 
æ SFP x MFP. 
Entry{A,Y¥}, modifies $5E,$5F,$62-$65,MFP,SFP,$AB. 

$E2B6: Multiply two byte integer tn $AD (low) and SAE (hich) with 
two byte integer in $64 (low) and Accumulator (high). Store 
product in X register (low) and Y register (high). 
Entry(A], Exit{Xelow byte product, ¥ «high byte product], 
modifies $65,$AD,SAE, puts $99 equal to zero. 

Or 

SDF4F: MFP = SFP OR MFP. Routine returns MFP = 0 if MFP = SEP 
æ 0, else routine retums MFP = 1. r+ 
Modifies MFP, puts $11 equal to zero. 

Power 

$EE97: MFP = SFPMPP, 
Entry(A = Z.F. = $9D}, modifies $T,$5E,$5F,$60-$66,358A-$8E, 
$92-$9C,MFP,$A3,SFP,$AB,$AE, SAD. 
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Subtract 


SE7AA: Determine $AB, subtract: MFP = SFP -- MEP. 


SE7A?: 


J. Nemee 


1. Purpose 


Modifies $92, MFP,SFP,$AB. 


Pull memory (> A,Y} in SFP, determine $AB, subt‘act: MFP = 
SFP - MFP. 


Entry(A,¥|, modifies $5F,$5F,$92,MFP,SFP,SAB. 


The Names routine can be used—during the evaluation of the & 
statement—to find the name, the status and the starting location of 
the sturage arca of simple variables, array elements and arrays. 


2. Name routine 


SDFEJ: At the start of the execution of $DFE3, the text pointer must 
point to the first character of the name. After the execution of 
SDFE3, the text pointer points to the first character behind the name 
cad the name and status locations are filled according to the table 


below. 
Name variable or array (el.) Status variable or array (el.) 
$81 $82 $11 $12 
heal pos pos 0 0 
Kring {S$} pos neg FF 0 
Bieger (%] neg neg 0 80 


For example, if a variable has the name AB, $81 and $82 will! 
contain the hex values 41 and 42 respectively, whereas if a 
variable has the name AB%, $81 and $82 will be loaded with 
the hex values C1 and C2. In the latter case, $12 is put equal 
to the hex value 80, to indicate that the variable is integer 
valued. 


Furthermore, Applesoft loads the pointer to the start of 
the storage area of the variable or the array in $9B [low] and 
$9C [high]. The pointer to the start of the storage area of the 
value of the variable or the array element is loaded in 
A=$83 (low) and Y =$84 {high}. If an array element is 
evaluated, the pointer to the first element of the array is 
stored in $94 (low) and $95 (high). 


In case one wishes to use matrix expressions in the & 
statement (for instance & A = A ~ B, where A and B are 
matrices], the hex value 40 must be stored in $14 before 
sxicuting SOFE3. Before returning tc BASIC, $14 has to be 
reset te zero again. 


Under the assumption hat no strings are used in the 
BASIC progran (which may leau to house cleaning 
activities), the tollowing locations may be modified by the 
execution of $OFE3. 


1. At the evaluation of simple variable names: 
$10, $11, $12, $81-$84, $94-$97, $9B, $9C, $B8, $B9. 

2. At the evaluation of array elements: $F, $10-$12, 
$81-$84, $94-$97, 39B, $9C, MFP, SAE, SAD, $B8, $B9. 
In addition, other locations may be modified, depending 
on the expressions in the subscripts. 


3. At the evaluation of {already dimensioned) array names 
which have to be interpreted as matrix names: $10, $11, 
$12, $81, $82, $9B, $9C, 3B8, $B9. 

{continued} 
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Apart from using names in the & 
Statement, one can also insert expres- 
sions. For instance, & SIN|1}) + SQR{3}. 
To evaluate such an expression, 
subrou.ine $DD67 <an be executed in 
the m.}. program. The result of the ex- 
pressiot is stored in the MFP. U the 
result has to be converted to an integer 
value, a JSR $E105 or a ISR S$E6ES in- 
struction can be used instead of the JSR 
$DD67 instruction. 


it might be possible that a wrong 
input to the mJ. program, or an error 
during the execution of the m.l. pro- 
gram, is detected. This will, for exam- 
ple, be the case if a matrix that is to be 
inverted turns out zo be not s square 
matrix. In that case, one may want 
to let Applesoft print an error 
message—indicating the kind of the 
error—wita tbe line number of the & 
statement that caused the etror. For 
this purpose, the routines listed under 
the heading Errors may be used. in the , 
case of the wrongly dimensioned i 
matrix, a JMP $E196 instruction, for i 
instance, displays the message "BAD 
SUBSCRIPT IN XX." After displaying 
the message Applesoft returns to 
BASIC command level. i 


Although there are more routines in 
table 1, it seems superfluous to discuss 
them here, since it will now be obvious 
how to use them. Instead, an example 
will be given to show how to integrate 
some of the routines in a matrix 
multipiication program. 


A Matrix Multiplication Program 


When written in BASIC, a matrix 
multiplication subroutine consists of 
the statements 


489 REM MATRIX MULTIPLICA. 
TION : C(R,P) = ARS) x 
ASP) 


500 
510 
520 
530 
540 


550 
560 
570 
580 
590 


FOR} 
FOA J 
LET D 
FOR K 
LET D 
B(K, |) 
NEXT K 

LET C(I} = D 
NEXTJ 
NEXT | 
RETURN 


3 a6 


oe = iaa 


oO 


ALI x 
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To execute this subroutine, the foliow- 
ing main program can be used: 


10 INPUT “DIMENSIONS 
MATRICES PRS? SPRS 

20 DIM AR,S)}BIS.P)OR,P) 
30 FOR! = 1 TOR 

40 FORJ = 1TOS 

50 LET AIJ = 14+ J 

60 NEXTJ 

70 NEXT i 

80 FORI = 1 TOS 
90 FORJ = 1TOP 
100 LET BOJ = 1 x 
110 NEXT J 
120 NEXT! 
130 GOSUB 500 
140 STOP 


J 


In the main program, the matrices A 
end B are dimensioned and initialized. 
Nex., tne mairix multiplication 
subroutine is called to put C equal to 
the product of A and B. If a m.l. pro- 
fam is written to multiply two 
matrices, the subroutine call at line 
130 can be replaced by 


IWéE&C=aAxEB 


Although the matrix names in the & 
statement can be chosen freely, we will 
use in the sequel the names C, A and B 
to denote the respective matrices. The 
dimensions of the matrices will be 
denoted by the same letters as in the 
BASIC program (i.e. P, R and S). 


The m.l. program can be split up 
into a main program and several 
subroutines. The main program per- 
forms the evaluation of the & state- 
ment. The first subroutine, called 
FNAME, takes care of the calculation 
of the pointers to the storage areas of 
the matrices C, A, and B in memory. 
The second subroutine, called 
WATMULT, is used for the actual 
matrix multiplication. Two other 
subroutines, ADD and ADDS, are call- 
w OY FNAME aan MATMULT to do 
some frequent!; occuring additions A 
discussion of the functions of the 
various foutines— which are listed in 
table 2— follows. 


1} The main program ($4000-$4022) 


The main program is written solely 
to control the multiplication of two 
matrices. Therefore it has to be re- 
placed by another main program if the 
number of matrix operations is extend- 
ed (with, for instance, add, subtract and 
inverse). The comment inserted in the 
listing shows how the program works. 


ee a Re ee ene ee ee ee 


Tabie 1: Applesoft Routines (continued) 


K., Normalize 


$E82E: Normalize MEP 
Exit(¥Y =O}, mcdifies MFP. 


L. Pack 


$EB72: Pack extension byte in MFP. 
Exit(X sentry, Y = entry), modifies MFP. 


Example < 
: JMP $5000 
: JSR $DFE3 
: JSR SLAF 
| JSA SDEBE 
. JSA $DFE3 


Jump to inultipiicaton program 


Find star:ing location of value of first variable 
Puli tirst variable into ihe MFP 

Check on comma in & statement 

Find starting location of value cf second 
variable 


` JSR | JSR $E97F Pull second variabie in the SFP, and 


: JSR $E982 Multic! MFP vatn SFP, store product in MFP 
: JSR SDESE Check on comma 

: JSA SDFE3 Find starting location of value of third variable 
: TAX Prepare entry store routine 

: JSR $EB28 Store product in third variable 


: ATS . Return to BASIC 





Table 2: Listings of Machine Language Programs 


A. The main program 


Purpose 


The evaluation of the & statement: &C a A x B, where C,A and B are 
matrices. 


Listing Comment 

$3FS : |MP $4000 Init jump location for & statement. 
$4000 : LDX #SFB Init location $08 for FNAME. 
$4002 : STX $08 
$4004 : JSR $4025 
$4007 : LDA #$D0 
$4009 : JSR SDECO 
$400C : JSR $4025 


] 
] 
$400F : LDA Sant 


Execute FNAME on first matrix {C). 
Check on '' =’ in & statement. 
Execute FNAME on second matrix (A). 
Check on ''x'' in & statement. 
Execute FNAME on third matrix {B). 


$4011 : JSR SDECO 
$4014 : JSR $4025 
$4017 : LDA $06 
$4019 : STA $71 
$401B : LDA $07 
$401D : STA $72 
S401F : JSR S40A5 Execute MATMULT. 
$4022 : RTS Return to BASIC. 


Restore column length of A {column 
length of C} in $71 and $72. 


. Subroutine FNAME 


Purpose 


Find name of array, check whether atray has two dimensions, zach less 
than 256. Store dimensions in $FC,X {second dimension) and SFD,X 
(first dimension}. Calculate column length of array (in bytes} and store 
it in $71 (low) and $72 (high). Calculate pointer to storage area of first 
element of second column ot artay, and store pointer in $0,X +2 and 
$1,X +2. FNAME can be called successively three times (or less}. Before 
the first call, the hex value F8 must be stored in location $08. A` the 
start of FNAME, the X register is loaded with the value in locatiol. 508. 
During the execution of FNAME the X register is incremented by two 
and stored in location $08 so that the contents of location $08 are in- 


cremented by two each time FNAME is called. kanini 
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2) Subroutine FNAME ($4025-$4083) 


Contrary to the man proyram, 
FNAME is constricted ia such a way 


that it can be used tor other matrix 
_ operations tou. The main purpose of 


e" ne Sp 42.04 


1 oR 


-- EF 


FNAME is to calculate the pinnter to 
the first element of the second column 
of the array being evaluated. The 
second column is taken because it is 
customary to use— when working with 
matrices—non zero values of the 
subssripte only, whereas Applesott 
leserves—when it encounter a DIM 
X(P,R} statement—P + | rows and R +l 
columns for the array because it allows 
ted subserpi As an example, sup 
pose that a DIM X23) instrucnon t9 
executed in a BASIC program. 
Applesoft stores the X-array column- 
wise {example 3, page 46}. 


“When multiplying the X matrix 


with another matrix, only the under- 


` scored elements have to be taken into 


Smee 


account. Since the first column con- 


tains no underscored elements, it can ` 


be skipped. 
Subroutine FNAME can be called 


_ successively (at most} three times. At 


y~ 


the first call, location $08 must contain 
the hex value FB, being the start of the 


: storage area, plus 4 of the matrix infor- 
< mation. Consulting the memory map 


Serre Tr 


of FNAME in table 2, it can be seen that 
the dimensions of the C array, P + 1 and 


` Rel, are stored in $F4 and SF5 and the 
' pointer to the first element of the 


second column of the C array |ie., the 


. pointer to C{(0,1} in SFA and SFB. Since 


eee 


Iecation $08 is automatically in- 
eremented by 2, each time FNAME is 

‘Hed, the dimensions of the next array 
2., the A amay} will—at the second 


_ teil of FNAME -—be stored in $F6 and 


$F7 and the pointer in $FC and $FD. 


: The information of the B array is stored 
- pSF8, $F9, SFE and $FF. 


ee 


<n Feria ee 


Apart fram the calculation of the 
pointer, FNAME also checks whether 
the array being evaluated has two 
dimensions, and whether the size of 
each dimension is less than 256. The 
lattes check is necessary because 
MATMULT can handle matrices with 
dimensions less than 255 only, which 
will be sufficient for almost all prac- 


- tical purposes. 


Finally, at each call of FNAME, the 
column length of the array being 
evaluated (which equals 5 times the 
pumber of column elements, since 


. Teals use S bytes of memory) is 
; Saleulated and stored in $71 and $72. 


Table 2: Lietings of Mactine Languagy Programs (continued) 


Memory Map of ENAME 


ees Column length (in bytes) of first array [C]. 

$08 Pointer to storage arca of array information. 

a Column length (in bytes) of third array {B}. 

$F4 Second dimension of first array (C}. 

$FS First dimension of first array [T]. 

$F6 | Idem for second array [A]. 

a ] idem tor third array (B) 

SFA Pointer to first element of second column of 

SFE first array {C}. 

Hae Idem for second array (A), 

ae Idem for third array (B). 

Listing FNAME Comment 

ae ; at . ety ] Put $14 equal to 40 for search of matrix name. 
$4029 : JSR $DFE3 Determine pointer to start storage area of artay. 
$402C : LDX $08 - X is loaded with pointer to storage area array 

information. 

nets i a di | Check whether array contains reals. 

$4032 : BEQ $4037 

$4034 : JMP $DD76 If not, display "TYPE MISMATCH". 

$4037 : STA $14 Put $14 back on zero. 

$4039 : LDY #504 

$403B : LDA #$02 

$403D : CMP ($9B),Y Compare number of dimensions of array with 2. 
$403F : BEQ $4044 i 
$4041 : JMP $E196 If not equal, display ‘‘BAD SUBSCRIPT”. 
$4044 : INY 

$4045 : LDA[$9B}, Y 

$4047 : BNE $4041 | Check whether dimensions of array are both less 
$4049 : INY than 256. If net, display '3AD SUBSCRIPT". 
$404A : LDA{S9B],Y| If yes, store second dimension in SFC,X and 

first dimension in $FD,X. 

$404C : STA $FC,X 

$404E : INX Note that at the first call of FNAME, the X 
$404F : INY register is loaded with F8. The dimensions of 
$4050 : CPY #$09 | the first array are thus stored in $F4 and SF5. 
$4052 : BNE $4045 

$4054 : TYA Accumulator contains 9 here. 

$4055 : CLC To obtain the pointer to the storage ares of the 
$4056 : ADC SOB first element in the array, ¥ is added to the 
$4058 : STA $00,X | painter in $9B and $9C The result is stored 
$405A : LDA $9C: in $00,X% flow) and $01,X (high). 

$405C : ADC #$00 | At the first call of FNAME X equals FA kere, so 
$405E : STA $U1,X J the pointer is stored in $FA and SFB. : 
$4060 : LDA $FB,X 

$4062 : LDY #$00 | Calculste the column length of the array by 
$4064 : STY $72 multiplying the size of the first dimension 
$4066 : ASL {which was stored in $FB,X) with 5. The colimn 
$4067 : ROL $72 length is stored in $71 and $72. 

$4069 : ASL 

$406A : ROL $72 

$406C : ADC $FB,X 

$406E : STA $71 

$4070 : BCC $4074 

$4072 : INC $72 
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$4074 : 
$4076: 
$4078 : 
S407A : 
$407C : 
$407E : 


$4081 


$4083 : 


CPX #$FA Isit the first call of FNAME! 


BNE $307E 

5TA $06 If yes, save column length of the array in $06 and 

LDY $72 $07. 

STY $07 l 

JSR $4087 Add column length to tke last calculated pointer 
: STX $08 to obtain the pointer to the storage area of the 

RTS first clement of the second column of the array. 


Store X in $08 for next calls of FNAME and 
return to main program. 


C. Subroutines ADD and ADD5 
Purpose ADD 


Add two byte integer in $71 (low) and $72 (high) to two byte integer in 
$00,X [uw] and $01,% {high}. Store result in $00,X (low) and $01,% 


(high) 


Entry|X}, Exiti/A= $01,X,X entry, ¥ = entry}. 


Purpose ADDS 


ADD 5 to two byte integer in $00, X {low} and $01,X (high}. Store result 
in $00,X (low) and $01,X {high}. 


Entry(X}, Exit{A= $00,X,X =entry,Y = entry). 


Listing ADD Listing ADDS 

$4085 : LDA $71 $4095 : CLC j 
i $4087 : CLC $4096 : LDA $00,X 
Foo‘ $4088 : ADC $00,X $4098 : ADC #$05 

$408A : STA $00,X $409A : STA $00,% 

$408C : LDA $72 $409C : BCC $40A0 

$408E : ADC $01,X $409 : INC $01,X 

$4090 : STA $01,X $40A0 : RTS 

$4092 : RTS 


D. Subroutine MATMULT 
Purpose 
Multiply matrix A{R,5} (dimensions in $F6 (S+1} and $F7 (r+1), 


witb matrix B(S,P} 


pointer in $FC and $FD} 


(dimensions in $F8 (P +1) and $F9 (S+ 1}, 
pointer in $FE and S$FF} 


and store result in 


matrix C{R,P} 


{dimensions in $F4 (P+ 1} and $FS (R+1]}, 
pointer in SFA and $FB) 


where P,R, and S each have to be less than 255. 


$06 
$97 
$17 | 
$18 
$19 
sia] 
$71 
$72 
$F4 
SF5 
$F6 


SFB 
$F9 
SEA 
SEB 
$FC 
$ED 
$FE} 
$FF J 
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Memory map of MaTMULT 


cr : Row counter for C. 
cs : Multiplication counter for row/ column multiplication. 


sp, : pointer to first element of current row of A. 
Pa pointer to first element of current column of B. 


k = S(R41}: Column length of A (in memory}. 


P + lLatentry. Used as column counter for C. 

R + 1 = number of elements per column of C {in memory) 
S+ 1: equals the number of multiplications necessary to 
multiply a row of A with a column of B. 


pc, : pointer to currert element of A. 
pce : pointer to current element of C. 
pa : pointer to first element of second column of A. 
pc, ; pointer to current element of B. 


(continued) 
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The column length of the C array, 
which equals the column length of the 
A array, is saved in locations $06 and 
$07, because the latter column | 

is needed later for MATMULT. 


3) Sudcoutine MATMULT 
($40A5-$4124) 


Before the matrices are multiplied, 
the dimensions ace checked to detet- 
mine whether thev satisty the condi- 
tions for muluplication. Next, the 
multiplication is carnied out as in- 
dicated by the fluw diagram in figure 1. 
The flow diagram shows that the pth 
column of C (p=2,..P+ 1} is obtained 
by multiplying the R rows of A jie., 
row 2,...R+ 1} each with the pth col- 
umn of B (p=2,...P + 1}. Note that at a 
row/column multiplication, the pro- 
duct of the first element of a row and 
the first element of a column is omitted 
since these elements have zero 
subscripts. The elements of the rows of 
A are separated by a distance of k 
(=S{R41}} bytes from each other in 
memory, so that each time a next row 
element of A is needed, k has to be add- 
ed to pca, After a row of A, not being 
the last row, has been multiplied with a 
column of B, the hpa pointer is in- 
cremented by 5 and pc, is put equal to 
bpa, so that pe, now points to the 
second element of the next row ot A. Lf 
a column of C, not being the last col- 
umn, has been filled, bpa is put equal 
to its starting value. That is, pa, and pg 
is put equal to pcp, whicb at that time 
points to the first element of the next 
column of B. 


The flow diagram further shows 
how “he multiplication of a row with a 
column is performed. The stack is used 
to store the sum of the products at in- 
ed so far, and each time a row elewtne 
is multiplied with a column element, 
the stack is pulled into the SFP and the 
newly-obtained product 1s added to it. 
The result is then pushed on the stack 
again. This process 1s continued until 
the row/column mulupiication t$ 
ready. The row/column product is then 
stored in memory. | 


Note that the subroutine address 
pushed on the stack at locations 
$40CE-$40D3 is the address minus one 
of the instruction following the jMP in- 
struction at location S40E5. If the 
MATMULET subroutine is relocated to 
another part of memory, this 
subroutine :ddress must be adjusted. 


a SE Se ee 


pe ee eT 


pou 
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Table 2: Listings of Machine Language Programe (continued) 


Liasing 


S4H0AS . 


SAA? - 
$4HHAY ; 


S40AR ` 


S40AE | 
$4080 - 
$4082.: 
$4084 ; 
$406: 
S40RB : 


S40HA : 
$40KC : 


$407 - 
SQu0KF : 
ha FAOS | 


$49CC 


$40CE: 
$40D0 : 


$40D1 


$40DE 
$40E0Q ; 
$40E2: 


$40E5 : 
$40E8 : 
$40EB : 


$40ED : 


$40F0 : 
$40F2: 
$40F5: 
s40F7 : 


SCOFS : 
SOFR : 


$40FD : 


200: 
42:02: 


$4105: 
$4107: 


$4109; 
$4108; 
$410E : 
$4119; 
$4112: 
$4114: 
$4115: 
$4118: 
$4114: 
$411: 
$4118; 
$4129: 
$4123: 


MATMULT 


LDA $t4 
CMI $EN 
N14) S4H0AE 
IM! $196 
LIDA $5 
CMP $i7 
RNE $IDAR 
LDA $ko 
CMP $E9 
BNF $40AB 
DEC $H4 
BNE $40BP 


RIS 
LIDA $ES ] 


: STA $06 
$40C3 : 
£40CE : 
asOC7 : 
$40C9: 
$30CA : 
: BM! $4100 


LDX #$03 ] 
LDA $FC,X 
STA §.7,X 
DEX 

BPL $40C5 


: LDA #$E7 
$40D3: 
$40D4: 
$40D? : 
$40D9 : 
$40DB : 


> LDA SFE | 


LDA #$40 
PHA 

PHA 

JSR SDEIO 
LDA $F8 
LDY $F9 
JSR SEAF9 


LDY SFF 
JSR $E97F 


JMP $DE47 
JSR SE7C] 
LDX #SF8 ] 
JSR $4085 
LDX #$FE ] 
JSR $4095 
DEC $07 
BNE $40CE 


LDY $FE 

JSR SEB2B 
LDX ¥SFA ] 
JE $4055 g 
Cec $06 
BEQ $4UBA 


LDX SFA l 


LDX +$)? ] 
ISR $4095 
STA SES 
LDA $18 
STA $F9 
LDA $19 | 
STA $FE 
LDA SIA 
STA SFF 
LDA $F6 ] 
STA $07 

JSR SES4E 


BEQ $40F0 _ 
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Comment 


Check dimensions for multiphcanon Ifan 
error is detected, display "BALD SUBSCRIPT", 


Pe P - 1: decrement column counter for. 
If equals zero, matrix multiplication Is 
ready 

Rerum to main program. 


ce Œ R + 1: init row counter for C. 


bpa = % 
Pa = CPi 


Always taker. 


Push subroutine address for routine $DE47 
on stack. 


Push MFP on stack. 
Load (pc,) in MFP. 


Load (pcp) in SFP and 

Mu:tiply MFP with SFP. Store product in 
MF?, 

Pul’ stack into SFP. 

Add MFP and SFP. Store sum in MEP. 


PC, œ pc, + k 


Pesa pea + § 


cs» cs - } 
If cs equals zero, row/column product is 
ready. 


Store MFP in (pcc}. 


Pec = pce +5 


cr = cr — 1 
is column of C filled? If yes, init hp,, p, and 
cr, 


hp, = hp, +5 
pc, = hp, 
PCa © Pa : restore column counter for B. 


cs = S + 1: init multiplication counter. 
Initialize MFP to zero. l 


Always taken. 
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Some Final Remarks 


Probably not all Apple owners will 
have Applesott in ROM. However, 
since the disk and tape versions of 
Applesoft which i have seen do sor dif- 
fer by more than a tew bytes itom the 
ROM version, i will be no big problem 
to convert the entries of the routines 
listed in table 1 to these versivus. The 
easiest way to do this 19 to find 
somcone whu hes Applesoft in ROM, 
80 that the diffezences can readily be 
traced back by comparing the versions 
with each other. In case nu ROM ver- 
sion is available, one can use the 
subroutine entry locations, which are 
found at the beginning of the Applesott 
Program. The sequence of the tirat 64 
subroutine entry locators’ cor- 
responds with the listing of the tokens 
in the Applesoft manual (#A2L0006 on 
pats 121). Next, the entry locations of 
the routines for SIGN to MIDS tollow. 
The rest of the entry locations’ are for 
+, ~, x, /,O, AND, OR, unary 
minus, NOT and comparison. Before 
each of the latter entries, a code, 
indicating the order of the operation, is 
inserted. 


Looking at table 3, where the 
entries of the routines for the ROM ver- 
sion are listed alphabetically, with the 
entries found in tape or disk versions of 
Applesoft, it will become apparent 
what differences there are. Aftez nat, 
the entry locations of the routine ¢ in 
table 1 can be converted accordingly. 


As a last point I wish to express my 
admiration for the ingenuity of the 
writers of Applesoft. During my study 
of Applesof:, I often searched tor bours 
for what was happening, ina seemingly 
endless sequence of {recursive} 
subroutines, which taught me a lot 
about m.l. programming. Apart trom a 
few errors that were made for instance, 
a zero byte forgotten between $E10) 
and $E102, so that the program 

10A% = -32768.00049 
PRINT A%e 


gives a surprising resule) ! came to the 
conclusion that Applesuft is a very 
good interpreter. 


I also wish to express my gratitude 
to Mr. F. Curvers of the Erasmus 
University in Rotterdam, for providing 
me with an excellent cross reference of 
Applesoft, without which my work 
would have been far more difficult. 


ee 
*Add one to the location found in the 
listing because the subroutines are 
executed via the RTS instruction. 
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Check dimensions Push sutroutine adress 


of matrices 





On stack 
Push FP? on stack 
Load (pc,) in ‘FP 


Load (pe, ) in SFP 
p cluleiply NF? vith S 


Pull stack in SFP- 


Add MFP and SFP 


ey RE ER 


(r 


Jute '!: F:o.t diagram of MAT,4 ULT. 


Remarks — See the memory map of MATMULT 
for the meaning of the symbols. 


- {x} mears: memory pointed xo by z. 
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Cornelis Bongers is an assistam: professor 
of statis:.cs at the Zrasmus University in 
Rotterdam. He uses his Apple U for solving 
statistical problems (for instance, 
likelihood maximatizetion}|. Another 
important field of application is finding 

' the solution of standardization problesas, 
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Good Articles 
and Good Photos 
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..« ‘dg increase in size—this 
ı .88u2 nas 96 pages, 18 more 
! than formerly—we can include 
Dieis AL. Wwidts. ol aO, wS ATF 
1 DRCN more S: Ecavs SL, 
„_ he articies we accept. 

MICRO is committed to 
covering all of the 6502 
systems. To do this well, we 
need a variety of articies on 
each system. We can always 
use more high-quality articies 
relating to AIM, SYM, KIM, 
Apple, Atari, PETIC8M, and Ohio 
Scientific systems. We are 
especially interested in good 
articles which apply to 6502 
systems in general. 


I Because we pian ‘co use 
more illustrations than formeily, 
we encourage authors to “think 
pictoriatiy” and to send us good 
tine drawings and black and © 
white photcs. 

‘Ne are aisc looking for black 
~—.nd whits phos which might 


: giang alow, ©." a def 
+ sapien. Pher.: of 6502 systacr.s 
‘ In ynusceal > sss 57 


. Sfofessionz. c. .0onsents 
«q SOUIE d8 espe ialiy wsicome. 
_ rhOtoz used i.ir sender ty ol 
. dvitelsa will be said for 
: separately. 


For details on how to submit 
manuscripts for possibie 
| publication, ask for MICRO 
ii Writer's Guide. Write or 
taiephone: 





Editorial Department 
MICRO 

®,3. Box 6502 
Ghe'msford, MA 01824 
#°7°98£-5515 
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X(2.3) 


which in exsence means: finding the set of 

sizes that minimizes the overall costs 

caused by the standardization of s product. 

Aa a hobby, he develope utility programs i 
for the Apple, such sa an assembler cross 

Plateada program and a dlsk-te-tape dump 

utility. 
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columr: 1 column 2 column 3 
aa ee ee Far 
T * fF fF + t t 7 
Xm, A0) X01) X21) X02) X22) X(0,3) X23 
j . 

pointer 

calculated 

by FNAME 





Teblo 3: Applesoft (ROM) Entry Locations | 





Engey 


$3F$ 
$E982 
$E7C1 
SETAA 
$EA69 


$SDF65 


SEBAF 
SOF55 
SEGES 
$FOSE 
$F105 
SE646 
$D66A 
or a4F 
$0896 
SEFEA 
$D995 
$E313 
$F331 
$DFD9 
$F769 
$0870 
$EF09 
$F280 
$E354 
$0766 
$E2DE 
SOBAO 
$0921 
$0935 
SFIS 
SF6E9 
SF3E2 


Dec Keyes Dec Key- Dec . 

tokes werd said token woi l meaty token py 
en F 

175 & $F3D8 144 HGA2 $F3BC 167 RECALL ' 

202 x $F286 163 HIMEM: $D9DC 178 REM : 

200 + $F232 142 HLIN $0849 174 RESTORE : 

201 - $FC58 151 HOME $F318 166 RESUME 

203 I $F6FE 147 HPLOT $0968 177 RETURN 

209 < SF7E7 150  HTAB S$E686 233 RIGHTS 

208 = $09C9 173 IF SEFAE 219 RNO 

207 > $FIDE 139 IN@  $F721 152 ROT= 

212. ABS $DBB2 132 INPUT $0912 172 RUN 

295 AMD  $EC23 211 . INT  $D8BO 183 SAVE 

230 ASC $F277 158 INVERSE $F727 153 SCALE =a 

197 AT $E65A 232 LEFT$ SODEF9 215  SCRN( 

225 ATN $E6DS 227 LEN $€B90 210 SGN 

140 CALL $DA46 170 LET  $F775 154 SHLOAD i 

a5" CHRE $D6A5 188 LIST  $EFF1 223 SIN j 

169 CLEAR $D8C9 182 LOAD $DB1E 195 SPC{ | 

159 SOLO T= $E941 220 LOC $£262 169 SPEED= ; 

187 CONT $F2A6 164 LOMEM: SEE8D 218 SOR i 

222 COS  $E691 234 MIDS — 199 STEP i 

131 DATA $0649 191 NEW  $D86E 179 STOP H 

184 DEF  $OCFe 130 NEXT S$F39F 128 STORE 

t23 DEL 452739 157 NORMAL $E3C5 226 STARS i 

134 DIM = $DE98 198 NOT $0818 192 TAR( 

148 DRAW S$F26F 156 NOTRACE S$FO3A 224 TAN 

128 END  $D9EC 180 ON $F399 137 - TEXT i 

221 EXP $F2CB 165 ONERR — 196 THEN Y tee, 

159 FLASH $OF4F 208 OR — £193 TO” ( 

19¢ FN $DFCO 218 POL  $F26D i55 TRACE j 

129 FOR  $E764 226 PEEK’ $A 213 USA 

214 FRE  $F225 141 PLOT $€707 229 VAL 

190 GET SE778 185 POKE $F241 143 VLIN 

176 GOSLƏ $D96B 161 POP $F256 162 VTAB 

171 GOTO $E2FF 217 POS $E784 18! WAIT 

136 GT $F1E5 138 PAS  $F76F 149 XDRAW 

146 ¥COLC = $D4DS 186 PRINT $EE97 204 A “3 

145 HGF COBE2 135 READ . ; 

| 
> >, 892 gu a 
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